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Abstract. In contrast to the situation with the 
geodynamo, no breakthrough has been made in the 
solar dynamo problem for decades. Since the ap- 
pearance of mean-field electrodynamics in the 1960's, 
the only really significant advance was in the field 
of flux tube theory and flux emergence calculations. 
These new results, together with helioseismic evi- 
dence, have led to the realization that the toroidal 
magnetic flux giving rise to activity phenomena must 
be stored and presumably generated below the con- 
vection zone proper, in what I will call the DOT 
(Dynamo-Overshoot-Tachoclyne) layer. The only 
segment of the problem we can claim to basically un- 
derstand is the transport of flux from this layer to the 
surface. On the other hand, as reliable models for the 
DOT layer do not exist we are clueless concerning the 
precise mechanisms responsible for toroidal/poloidal 
flux conversion and for characteristic migration pat- 
terns (extended butterfly diagram) and periodicities. 
Even the most basic result of mean-field theory, the 
identification of the butterfly diagram with an et- 
ui dynamo wave, has been questioned. This review 
therefore will necessarily ask more questions than 
give answers. Some of these key questions are 

- Structure of the DOT layer 

- a-quenching and distributed dynamo 

- High-latitude migration patterns and their in- 
terpretation 

- The ultimate fate of emerged flux 



1. INTRODUCTION 

The turn of the millennium invites us to look back 
and draw balances in all fields of human activity. Yet 
in solar dynamo theory we also have an added incen- 
tive to make such an assessment. In the theory of 
the geodynamo a significant breakthrough has been 



achieved in the past few years ( Glatzmaier 


fe Roberts 


1995 




Kuang & Bloxham 


1997 


, plson et 


al 




1999), 



One cannot but wonder if a similar breakthrough is 
within reach in the case of the solar dynamo. Un- 
fortunately, as it will turn out from this review, the 
prospects are rather bleak, at least on a short term. 

As for such a comparative assessment one needs a 
wider historical outlook this review will not be re- 
stricted to the developments that have taken place 



leading to a surge of renewed activity in the field. 

In: The Solar Cycle and Terrestrial Climate, 
ESA Publ. SP-463, p. 3-14 (2000) 



since the reviews of Weiss (1994) and Schmitt (1993). 
(Such developments were mostly limited to advances 
in the study of interface dynamos, cf. Sect. 3 be- 
low.) A wider historical overview, starting with the 
dawn of mean-field theory in the 1950s and 60s will 
thus be given in Section 2 below. Given the finite 
amount of space available, I will compensate for this 
wider temporal scope of the review by restricting at- 
tention strictly to the problem of the origin of the 
solar cycle, i.e. of the 22 year periodic variation of so- 
lar activity, and associated migration patterns (but- 
terfly diagram). Solar activity variations on both 
shorter and longer timescales are ignored, as are 
the solar-type magnetic cycles in other stars, non- 
axisymmctric phenomena such as active longitudes, 
and the problem of the long-time phase coherence of 
the cycle. This restriction is imposed out of necessity 
only, and in no way does it imply that these effects 
do not yield important clues even to the origin of 
the 22 year cycle itself. Clearly, a critical test of any 
theory of the solar cycle is whether it can be readily 
extended to predict these other phenomena as well. 

After the historical overview, Section 3 will at- 
tempt to cut some order in the dazzling multitude 
of solar dynamo models by introducing a classifi- 
cation scheme. Three main model families can be 
clearly discerned: overshoot dynamos, interface dy- 
namos and flux transport models, circulation-driven 
"conveyor belt" models being the most important 
subgroup of the latter class. Finally, Section 4 calls 
attention to some key areas where more intensive 
theoretical or observational efforts could lead to sig- 
nificant advance. 

But first of all we should state clearly what are the 
basic observational facts to be interpreted by a so- 
lar cycle model. Once we apply our aforementioned 
restriction excluding long-term variations, stellar ac- 
tivity etc., the remaining list is quite short. 
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The 11/22 year cycle period. Beside re- 
producing the value of the period, the crude 
agreement of this value with the timescale of 
pole-equator diffusion in the convective zone 
also asks for an explanation. While such an 
order-of-magnitude equality can certainly be 
coincidental (cf. the coincidence of the solar ro- 
tation period with the convective turnover time 
in the deep convective zone), a natural explana- 
tion for it would clearly make any cycle model 
more attractive. 
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Figure 1: Extended butterfly diagram of solar activity: 
time-latitude distribution of sunspot groups (low-latitude 
branches) and polar faculae (high-latitude branches). Af- 
ter Makarov & Sivaraman (1989) 

- The characteristic migration pattern 
(extended butterfly diagram) . Our knowl- 
edge of latitudinal migration patterns of mag- 
netic flux is summarized in the extended but- 
terfly diagram of Figure [l]. The tracers shown 
here track partly toroidal and partly poloidal 
fields.^ While azimuthal field lines by defi- 



nition cannot cross the surface, the observed 
properties of large-scale solar active regions^ 
strongly suggest that they are tracers of a sub- 
surface toroidal field, locally bulging out into 
the atmosphere. In this sense, photosphcric 
magnetometry can give us information about 
the migration patterns of the toroidal field 
component as well. At low 35°) latitudes 



2 Note that, owing to our assumption of axial symmetry, 
"toroidal" is now synonymous with "zonal" or "azimuthal", 
denoting the </>-component of a vector field in spherical co- 
ordinates, while "poloidal" is synonymous with "meridional", 
denoting the remaining components. This relieves us from 
giving a more generic definition of these terms. 

3 To be specific: their preferential East-West orientation 
following Hale's polarity rules. (The western, or leading 
magnetic polarity is identical on the same hemisphere and 
within the same cycle, and alternates between hemispheres 
and cycles) 



both the poloidal and the toroidal field compo- 
nents migrate equatorward. At high latitudes, 
poloidal fields show a marked poleward migra- 
tion, as indicated also by the migration pattern 
of a number of tracers such as quiescent promi- 
nences or the coronal green line. The migration 
pattern of high-latitude toroidal fields is less 
clear — a point we will return to in Section 4.2 
below. 

The confinement of strong activity (large 
active regions) to low heliographic lati- 
tudes |$| < 35°. 

The phase dilemma(s). In its original sense 
( [5 fix 1976| ) the phase dilemma consists of the 
fact that at low latitudes the radial field (de- 
rived by azimuthal averaging of the magne- 
tograms) is in an approximately 7r phase lag 
compared to the toroidal field at the same lat- 
itude. Another phase lag to be explained is 
the 7r/2 lag between the two branches of the 
butterfly diagram, i.e. that the polar field re- 
versal occurs slightly after the sunspot maxi- 
mum. Finally, the phase of torsional oscilla- 
tions (cycle-related periodic oscillations of the 
rotational velocity in migrating belts) relative 
to the toroidal field is a third quantity con- 
straining theories of the cycle. 

2. HISTORY 
2.1. Convection zone dynamos 



It all started with Parker's ( 1955| ) classic paper that 
set down the foundations for solar aft dynamo the- 
ory. In its trace, mean fiel d electrodynamics was 
developed during the 1 960's ( Bteenbeck et al. 19(iC| . 
Krause fc Radler] |l980| ). To give a reminder of the 



basics, described in so many other reviews (e.g. Cowl- 
ing 1981 , Belvedere 1985 ) the induction equation in 
a turbulent conductive medium reads 

<9 t B = Vx(UxB + £)-Vx ryV x B. (1) 

where B is the mean magnetic field, U is the large- 
scale flow velocity, dt denotes time derivative, rj is 
the magnetic diffusivity, and £ is the turbulent elec- 
tromotive force arising as a result of the interaction 
between the turbulent velocity field v and the turbu- 
lent magnetic field. This latter is in turn the result 
of the action of v on the mean field B, so £ is a 
functional of v and B. Assuming scale separation 
I <C Hb where Hb is the length scale of the mean 
field and I is the scale of turbulence, £ can be ex- 
panded in the derivatives of B. For homogeneous 
and isotropic turbulence this yields 



£ = aB-/3VxB 



(2) 
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where a and (3 are now functionals of v only. Substi- 
tuting (||) into (|l|) we see that the role of (3 is formally 
identical to that of rj. For this reason (3 is called tur- 
bulent magnetic diffusivity, and elementary consid- 
erations or even dimensional analysis yield (3 ~ Iv. 
As for turbulence the Reynolds number /3Jrj ^> 1, in 
practice ij can be omitted in equation (pi). In con- 
trast, the pseudoscalar a gives rise to a qualitatively 
new effect, the a-effect. 

In the axisymmetric case considered here, using 
spherical coordinates 8, (f>, r, B can be split as 

B = Be + V x (Ae ) 

where is the azimuthal unit vector, B is the 
toroidal field component, and A is the (toroidal) vec- 
tor potential of the poloidal field. We further assume 
that U is a pure rotation 

U = r sin Oile^ 

and introduce the shear vector 

n = r sin 9 VO 

In the limit rsinO/Hg — > oo the form of the vector 
operators simplifies to their form for the local Carte- 
sian frame e x > — eg, e y i — e^, e z i — e r . Now we 
introduce a new frame xyz by rotating x'y' z' around 
e^ with an angle — ir/2 < A6> < tt/2 so that Cl x — 0. 
With these assumptions and notations the poloidal 
and toroidal parts of equation (Q) read 



d t A = aB + 13V 2 A 



dtB = fl z d x A 



aV 2 A 



(3V 2 B, 



(3) 
(4) 



known as the classic dynamo equations. 

It is clear from (^-^) that the role of the pseu- 
doscalar a is to turn the poloidal and toroidal field 
components into each other which implies some kind 
of helical motion. The classic candidate for this, sug- 



gested by Parker (1955) is the passive advection of 
fields by helical convective motions. Later, alterna- 
tive mechanisms for an a-effect were also proposed, 
based on a dynamic interaction of field and motions 
(see Section 4.3 below). A general property of these 
mechanisms is that a turns out to be positive in 
the bulk of the solar convective zone in the north- 
ern hemisphere, while it tends to be negative in the 
stably stratified layer below. (Being a pseudoscalar, 
a changes sign between hemispheres.) 

The shear fl, associated with differential rotation, 
in turn, winds up the poloidal field into a toroidal 
component. Without the a and Q terms we would 
be left with a diffusively decaying field, so at least 
one of these terms is necessary for dynamo action in 
both equations. Depending on which, if any, of the 
dynamo terms in (0) is discarded, we distinguish a 2 , 



afl and a 2 fl dynamos, a 2 dynamos can be shown 
to give rise to non-oscillatory behaviour and toroidal 
and poloidal field amplitudes of the same order of 
magnitude which does not agree with the properties 
of the solar dynamo. Thus, in what follows we will 
concentrate on ail dynamos, neglecting the second 
term on the r.h.s. of (|]). (Note, however, that under 
more general conditions than those considered here, 
oscillatory a 2 dynamos can also be constructed , as 
pointed out recently by Schubert fc Zhang 2000 .) 

Assuming that beside a and f3, fl z can also be 
regarded constant, the system (^-|J) is homogenous 
and linear, admitting wavelike solutions of the form 

B = B Q exp [i(ut - kx)] (5) 

A = A exp [i(ut - kx + 5)] (6) 

where lu = ujr + iluj is a complex frequency while 
all other variables are real. Wr, Aq and Bq can be 
taken to be non-negative without loss of generality. 
Introducing the (signed) Reynolds numbers 



a/ /3k R n = n z /pk 2 



(7) 



as well as the dynamo number and the nondimen- 
sional frequency 

D = R a R n u = uj/[3k 2 (8) 

and substituting the Ansatz (|]-§]) into (||-[|) we find 

{l + iu)A = kr 1 R a B () e- lS (9) 

(1 + iu>)B = -ik x R n A e u (10) 
The product of these latter equations is 

(1 + i£o) 2 = -iDk x /k (11) 

As only k x appears in (|l^), no unstable modes 
(self-excited field or "dynamo waves") exist with 
k x = 0. Remembering the way we oriented our x- 
axis this implies that dynamo waves propagate along 
isorotational surfaces. But in which direction? The 
imaginary part of ([ll]) reads 



2w_r(1 - Qi) = -Dk x /k 



(12) 



As for unstable modes (self-excited field or "dynamo 
waves") uji < 0, (|l2] ) yields the important result 



k x D < 0, 



(13) 



known as the Parker Yoshimura sign rule. Thus, 
e.g. in the northern hemisphere equatorward propa- 
gation (k x > 0) implies D < 0. With a positive a, as 
is the case in the bulk of the convective zone, this im- 
plies d r fl < 0, an outward decreasing rotational rate. 
This was indeed the general expectation for the solar 
internal rotational law in the 1960's and 1970's. 
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The solution of (|llj) is 

u> = i±{l-i){-Dk x /2k) 1/2 



(14) 



For unstable modes obviously the plus sign applies. 
So the growth rate is 



= {-Dk x /2k) 1/2 - 1 



(15) 



Unstable modes thus exist when \D\ > 2 in nondi- 
mensional units. As according to equation (||) \D\ 
decreases with k, it is the lowest k modes, with a 
scale comparable to the solar radius Rq that have 
the highest growth rate and will dominate the so- 
lution. (Note that this implies that our formalism, 
derived for the limit kr sin 9 — > oo, is strictly speak- 
ing invalid for these modes — nevertheless it may still 
be used for general guidance.) 

Note that when Cuj = 0, lor = 1 follows from (Q): 
the period of the critical mode is thus just the 
diffusive timescale corresponding to Rq. Esti- 
mating a and tt z on the basis of helical convection 
and the observed differential rotation, D proves to be 
order of unity for a convection zone dynamo, show- 
ing that the dynamo is indeed approximately critical, 
and thus naturally explaining the agreement of the 
cycle period with the diffusive timescale. 

Finally, let us note that it is straightforward to 
work out from the above formulae that for an equa- 
torward propagating wave, the phase of the radial 
field component relative to the toroidal field is — Stt/A 
if a > in the northern hemisphere and 7r/4 other- 



wis^JStix 19761 



Taken altogether, the above considerations showed 
that for the expected positive a-effect in the solar 
convective zone, assuming an inwards increasing ro- 
tational rate, one can correctly reproduce the cycle 
period, the equatorward branches of the butterfly di- 
agram as a dynamo wave, and the low-latitude phase 
relationship. (The high-latitude poleward branch 
could obviously be reproduced by assuming d r Q > 
there, though this line was not pursued, relying on 
the Babcock-Leighton approach instead.) All this 
gave the impression that, missing details apart, the 
basic mechanism of the solar dynamo is well under- 
stood. 

2.2. Crisis 

The first warning signs that something is amiss be- 
gan to appear towards 1970, with the realization that 
most of the magnetic flux in the solar photosphere, 
and presumably below, is present in a strongly in- 
termittent form, concentrated into strong flux tubes 



(Weiss 1964, gheeley 1966, Stcnflo 1973, Howard fc 



Stetiflo 1972). Flux tube theory developed in the 
1970's and it became clear that for thicker flux bun- 
dles, owing to their lower surface/volume ratio, vol- 



ume forces such as buoyancy, curvature and Corio- 
lis forces dominate over the drag of the surrounding 
plasma flows acting on the surface. These tubes can 
then move largely independently of the surrounding 
flows, invalidating simple one-fluid descriptions like 
mean field theory. Thus, solar magnetic fields can 
be divided into two components: passive fields, con- 
sisting of thin flux fibrils, that move passively with 
the flow owing to the drag and are the subject of 
mean field theory; and active fields, consisting of 
thick flux bundles moving under the action of vol- 
ume forces. And the characteristics of large active 
regions strongly suggested that they are essentially 
(fragmented) loops formed on toroidal flux bundles 
of 10 22 Mx which clearly fall in the active category, 
outside the jurisdiction of mean field models. 



What is more, Parkei ( 1975 ) called attention to 
the fact that such flux bundles cannot be stored in 
the convective zone for a time scale comparable to 
the cycle period, being subject to buoyant instabil- 
ities that can rapidly remove the whole tube from 
the zone. The only place to store these tubes is near 
the bottom of the convective zone, especially in the 
stably stratified but still turbulent lower overshoot 
layer below it. Toroidal flux tubes lying here may 
still develop finite-wavelength buoyant instabilities 
that may give rise to loops erupting through the con- 
vective zone into the atmosphere, producing active 
regions. This scenario has gained firm foundations 
with the first nonlinear calculation of the emergence 



of such loops through the convective zone ( Moreno 
Inscrtis) 198E), and such flux emergence models have 



by now evolved into an independent chapter of the 
global dynamo problem. In the present review we 
will not deal with this topic in detail (see the review 



by Moreno- Insertis 1994 ), even though flux emer- 
gence models are the only real "success story" of dy- 
namo theory since the 1960's. While many details 
are still unclear, by now these models can reproduce 
sunspot proper motions and active region flux distri- 
butions to a quite convincing detail. A very robust 
main conclusion from the models, of great impor- 
tance for the global dynamo, is that in order to re- 
produce the observed characteristics of active regions 
the toroidal flux tubes must have a field stregth of 
about 10 5 G — an order of magnitude higher than 
the turbulent equipartition field in the deep convec- 
tive zone. Explaining the origin of such strong fields 
is a major challenge for dynamo theory. 

Guided by these realizations, in the 1980's the first 
attempts were made to construct dynamos operating 
in the lower overshoot layer. The unknown profiles 
of a and the differential rotation, however, greatly 
impeded progress, allowing a far too wide parameter 
space to play with. Therefore, attempts were made 
to numerically simulate the whole convective zone, 
with a consistent picture of differential rotation, heli- 



4 



cal convection, and dynamo flGlatzmaici 1985 ). Nev- 
ertheless, the results (poleward propagating dynamo 
waves) were at odds with the observations, and when 
finally even the predicted differential rotation profile 
(constant on cylinders) was proven wrong by helio- 
seismic measurements (constant on cones), the sim- 
ulational approach was abandoned (cf. the remarks 
in Section 5). 

On the other hand, the helioseismic determination 
of internal differential rotation gave new impetus to 
mean field dynamo theory. Those inversions clearly 
showed that most of the shear (the fi-effect) is con- 
centrated in a thin layer near the bottom of the con- 
vective zone, known as the tachocline. This was seen 
as further evidence that a thin layer situated about 
20QQ0Qkrn below the solar surface is nf \re\ 



Mil 
uJc 



°y impor- 



tance for the working nf the solar dynamo Depenrl- 



ing on the physical viewpoint we study it from, this 
layer is alternatively called dynamo layer, overshoot 
layer or tachocline. In the present review for sim- 
plicity I will refer to it as the DOT (Dynamo - 
Overshoot — Tachocline) layer. 



Table 1: A classification scheme for solar dynamo mod- 
els 
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3. MAIN FAMILIES OF MODELS 

Solar mean field dynamo theory in the past decades 
gave rise to a bewildering variety of models. Never- 
theless, all recent models (i.e. those using the he- 
lioseismic rotation law) can be classified into just 
four main types according to a plausible classifica- 
tion scheme (Table [I] and Figure ||). One classifi- 
cation parameter here divides the models according 
to whether they still interpret the butterfly diagram 
as a dynamo wave just like the orthodox convection 
zone dynamos did, or if they substitute that with 
some flux transport mechanism (meridional circula- 
tion or pumping). The other parameter in turn di- 
vides the models according to whether the a- and f2- 
effects are cospatial or "distributed", i.e. they take 
place in different (adjacent or very distant) parts of 



the model volume. The resulting four model types 
are as follows. 

3.1. Cospatial wave models: OL 
dynamos 

Widely known as "overshoot layer" or OL dynamos, 
these are perhaps the most conservative models that 
simply replant the concepts of the convective zone 
dynamos of the 1960's and 70' into the DOT layer. 
As the helioseismic inversions show <9 r fi > at low 
latitudes, these models need to assume a < to get 
the right migration directions. This assumption is 
rather plausible in the DOT layer for several different 
physical mechanisms for a. The state-of-the-art in 



this approach is r epres ented by the model of Rudiger 



Brandenburg ( 1995 ) 



Successes: 

- The butterfly diagram comes out right. The 
Parker- Yoshimura rule leads to polar and equa- 
torial branches separated at the corotational 
latitude, as observed. 

Difficulties: 

- The low-latitude phase dilemma: as now a < 
0, the radial field is found to be nearly in phase 
with the toroidal field, instead of being in an- 
tiphase as observed. 

- The cycle period tends to be too short for 
thin layer models. This is basically because 
the same amount of differential rotation is now 
concentrated in a much thinner layer, leading 
to much stronger shear and much higher dy- 
namo numbers. This may be compensated by 
reducing a, e.g. by arguing that, as nonlin- 
ear effects act via a quenching of the a and 

mechanisms, it is natural to expect that af- 
ter saturation the dynamo will be effectively 
critical. Yet the degree to which the nonlin- 
earity can increase the cycle period depends a 
lot on the assumptions made, and in general 
it does not seem sufficient (Rudiger 1994). At 
any rate, certain "tricks" , such as using an a 2 VI 



dynamo (Oilman et al. 1989) or introducing an 
intcrmittence factor (Riidigcr & Brandenburg 



1995), can save the models, but then the crude 
coincidence of the period with the lateral dif- 
fusive timescale is coincidental. 

For an equatorial confinement of strong fields 
the a-effect needs to be arbitrarily confined 
to low latitudes. Nevertheless, such a- 
distributions are indeed found in some calcu- 



lations (cf. gchmitt 1993). 



3.2. Distributed wave models: IF 
dynamos 




Meridional circulation 



Pumping 

S Buoyant emergence 

Figure 2: Schematic illustration of the main features of main solar dynamo families. The sites of a and ft effects 
and toroidal/poloidal field transport mechanisms and directions are shown in a meridional quadrate. The internal layers 
marked schematically are, from the inside outwards, the radiative, DOT and convective layers. OL: Overshoot layer 
dynamos; IF: Interface dynamos; CP: Cospatial pumping models; BL: Babcock-Leighton type models; 
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Parker ( 1995 ) suggested a dynamo where the diffu- 
sivity discontinuously varies by many orders of mag- 
nitude across a surface. The a-effect operates on 
the high-diffusivity side of the interface, while the f2- 
effect (shear) is limited to the low-diffusivity side. He 
showed that under these conditions a dynamo wave 
can be excited, obeying the Parker- Yoshimura sign 
rule. The attractive feature of this model is that the 
toroidal field generated on the low diffusivity side can 
be made arbitrarily strong by reducing the value of 
magnetic diffusivity there. Thus, the origin of 10 5 G 
fields could be explained. 

Parker's analytic, plane parallel model has been 
extended to more realistic situations and incorpo- 
rated in full solar dynamo models in a number of pa- 



pers (Charbonneau 
Th 



MacGregor 1997, Markiel fc 



1999). Unfortunately, the results are some- 



what contradictory owing to numerical problems re- 
lated to modelling the discontinuity. In this context 
we may perhaps note that the main physical differ- 
ence of these IF models compared to the OL mod- 
els is the spatial separation of a and Q. The intro- 
duction of a discontinuity between them is an added 
feature that can simplify the analytic treatment but 
at the same time complicate the numerical calcula- 
tions. A model where the diffusivity is continuously 
distributed, albeit with a sharp gradient, with a and 
n concentrated on the two sides of this gradient, may 
well be worth considering. Beside being more real- 
istic, such a model might also avoid the numerical 
problems mentioned. 
Successes: 

- Strong toroidal fields can be readily explained. 
Difficulties: 

- At the present stage, numerical difficulties pre- 
vail. An evaluation of the physical performance 
of the models will only be possible after these 
are resolved. 



3.3. Cospatial transport models: CP 
dynamos 

In an inhomogeneous medium, the proper expression 
of £ is more general than equation (||), the scalars 
a and [3 being substituted by tensorial expressions. 
A tensorial d-term in (^|) can then alternatively be 
written as 

dB = aB + d<sB — 7 x B 

where a — Tr a, the last term is the vectorial product 
equivalent of the action of the antisymmetric part of 
d on B, and as is the symmetric and traceless part of 
d. Substituting (^) into ([j]) we see that the role of —7 
is analoguous to that of U, i.e. it formally describes 
an advection of the magnetic field. This effect is 
called the (normal) pumping of the field along the 



inhomogeneity. The ag-term can be shown to give 
rise to a similar effect with the difference that the 
sign of the transport depends on the orientation of 
the field component perpendicula r to the pump ing 
direction (anomalous pumping; see Petrovay 1994 for 
a detailed discussion). 

It must be stressed that the descriptions using a 
tensorial alpha and those using pumping effects are 
formally equivalent. The advantage of the pumping 
formalism is that it helps one to get a physical "feel- 
ing" of the processes at work. 

Depending on the particular inhomogeneity asso- 
ciated with the pumping, one can speak of density 

pu mping, turbulent pumping e tc. It was suggested 

by Krivodubskij fc Kichatinov ( 1991 ) that, instead 



of appealing to a dynamo wave, the field migration 
patterns can also be interpreted by density pumping, 
directed towards the rotational axis for poloidal fields 
and away from it for toroidal fields. This pump- 
ing is supposed to operate throughout the convective 
zone, instead of being confined to the DOT layer. A 
more detailed model along these lines has recently 
been constructed by Kitchatinov et al. ( 1999] ). While 
the model reproduces well the phase relation for the 
torsional oscillations, it does not even address ques- 
tions such as the origin of the deep-seated toroidal 
field. Thus, while the basic concept is interesting, at 
the present stage this "cospatial pumping" approach 
cannot be regarded as a serious aspirant for the ex- 
planation of the global dynamo. 

3.4. Distributed transport models: BL 
dynamos 

These are generally known as "Babcock-Leighton- 
type" models (hence the BL code). While they in- 



deed grew out of the semiempirical approach of Bab 



cock (1961) and [Lcighton ( 1964| ) to the solar cycle 



they did get a lot more radical in the 1990's. In- 
deed, in the era of convective zone dynamos, in the 
1960's and 70's, the Babcock-Leighton approach was 
not seen as necessarily conflicting with the dynamo 
wave theory of the cycle. With a dynamo operating 
throughout the convective zone, the a-effect caused 
by the inclination of active region axes relative to E- 
W (which is the cornerstone of this approach) could 
be considered as helical convection caught in the act; 
and in h i s mat hematical formulation of the model 
Leightorj ( |l964| ) explicitly assumed that the toroidal 
field migration is the result of a dynamo wave. It 
was only the poleward migration of the poloidal field 
that was interpreted in terms of a lateral transport. 

In subsequent versions of the model, meridional 
circulation plays the main role in transporting the 
poloidal fields to the poles near the surface. And 
from here it was just one step to close the circle and 
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assume that the deep return flow of meridional cir- 
culation is responsible for the equator ward drift of 
the toroidal field. This step, made by Wang et al 



( 1991 ) was radical indeed: it represented the first 
break with the canonical interpretation of the but- 
terfly diagram as a dynamo wave, generally accepted 
for three decades. 

The model essentially works like a conveyor belt: 
the poleward meridional circulation near the surface 
transports the poloidal fields towards the poles at 
high latitudes, giving rise to the poleward branch of 
the butterfly diagram. At the poles, the fields are 
advected down to the bottom of the convective zone 
where the shear converts them into toroidal fields 
that get amplified while advected towards the equa- 
tor. Once these are strong enough, they are sup- 
posed to form buoyantly emerging loops at low lat- 
itudes that give rise to active regions, the Coriolis 
force lending an inclination to the loop planes, i.e. 
introducing an a-effect, and thus regenerating the 
pol oidal field. 



f ko re re c ent versions of the model ( C houdhuri| 



et |l995 , Dikpati fc Charbonncau 1999 ) have de- 
veloped it into an internally consistent modelling ap- 
proach with the ambition of yielding a complete de- 
scription of the solar cycle. As such, this family of 
models is at present the most serious competitor of 
the OL class. 
Successes: 

- The low-latitude confinement of strong activity 
comes out rather naturally. (The toroidal field 
is amplified by shear as it is advected equator- 
ward.) 

- A tolerable reproduction of the extended but- 
terfly diagram (although the latitude where the 
two branches part tends to be too high). Note, 
however, that the original (one-dimensional) 
Babcock-Leighton models gave a much closer 
fit to the observed migration patterns, which 
was their main asset. In the extension to two 
dimensions, BL dynamos had to sacrifice this 
achievement. 

Difficulties: 

- The agreement of the corotational latitude 
with the latitude where the butterfly branches 
diverge must be considered coincidental (if re- 
produced at all). 

- The a-effect is confined to the surface. This is 
probably unrealistic even for flux tube alpha, 
cf. the discussion in Section 4.3 below. 

- The model, as all models relying on a flux tube 
alpha, is not self-exciting, as this a-effect works 
for strong fields only. It needs another dynamo 
mechanism to "kick it in" . 



A serious problem is that the approach only 
works with an unrealistically low value for the 
turbulent diffusivity in the convective zone, 
(3 ~ 10km 2 /s. (This is needed to keep the 
two parts of the "conveyor belt" separated. 
In fact the main difference between the BL 
and IF models is that the separated seats of 
a and f2 communicate by diffusion in the IF 
models.) There is no justification for such 
a low value. On physical grounds, (3 ~ 
Q.blv ~ 500km 2 /s is expected — indeed, even 
the one-dimensional Babcock-Leighton models 
predicted 600km 2 /s. According to the contro- 
versial proposal of Vainshtein & Rosner (1991), 



on the other hand, the diffusivity should have 
been quenched by many orders of magnitude to 
values far lower than those used in the BL mod- 
els. (Note that by now this question has been 
settled: diffusivity suffers no strong quenching 
in realistic, intermittent fields in three dimen- 



sions, Gruzinov & Diamond 1994, Petrovay & 



Zsargo 1998 
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Figure 3: Steamy window: with realistic turbulent trans- 
port parameters, an arbitrary time- dependent poloidal 
field pattern imposed at the bottom of the convective zone 
is reflected, somewhat blurred, at the surface. After 
Petrovay & Szakdly (1999) 

One may wonder if the effect of the strong dif- 
fusive link can be reduced by some kind of se- 
lective field pumping that would keep the two 
field components apart 



Petrovay & Szakaly 



(1999) investigated this possibility in a model 



of poloidal field transport, incorporating all 
known transport effects (circulation, diffusion, 
turbulent pumping, density pumping). It turns 
out that diffusion is dominant, forging such a 
strong link between the DOT region and the 
surface that any migration pattern imposed at 
the bottom will be reflected at the top (Fig. ||). 
This seems to represent a serious difficulty for 
BL models. 

4. KEY ISSUES 

4.1. Structure of the DOT region 

It is obvious from the above model descriptions that 
a thorough understanding of the structure and dy- 
namics of the DOT layer is crucial for advance in so- 
lar dynamo theory. Unfortunately, at present we do 
not have a reliable model for this layer. As a compli- 
cated interplay of convection, rotation and magnetic 
fields is expected there, as a first step one would ex- 
pect to develop models for just one aspect of the 
problem: a model for overshoot with no rotation 
and magnetic fields; a non-magnetic model for the 
tachocline with no overshoot; etc. Nevertheless, we 
lack even such simplistic models, and surprisingly lit- 
tle effort is made to construct them. 

In the case of overshoot, the way to construct 
such a model is in principle well known: it consists 
in solving a conveniently truncated hierarchy of the 
Reynolds stress equations. Marik fc Petrovay ( [2000 ) 
present preliminary results from a numerical study 
of this problem. 

Helioseismic inversions show that the thickness 
of the tachocline is about 0.1 Rq. This is puz- 
zling, as even in the absence of turbulent diffu- 
sion, the Eddington-Sweet circulation should have 
mixed these layers enough in 4.6 • 10 9 years to ex- 



tend the tachocline to much greater depths (Spiegel 
& !5ahB 1992p . It is most often assumed that the 
solution to this "thin tachocline problem" resides 
in anisotropic turbulence: an extremely strong hor- 
izontal anisotropy would indeed effectively smooth 
out horizontal gradients without contributing much 
to the vertical transport. But in view of strong 
nonlinear mode coupling in turbul ence, such an ex- 
tre me an is otrop y seems dubious. Forgacs-Dajka fc 
Pel rovay ( [2000 ) explore the alternative possibility 
that the necessary horizontal momentum transfer 
is due to an inverse A-effect: however, the neces- 
sary amplitude of A proves even more unrealistically 
large. Magnetic decoupling of the envelope from the 
radiative interior is a third explanation put forward 
by jRudiger fc Kitchatinov] (|1997|). 



models, one needs an issue on which they yield con- 
flicting predictions that are relatively easy to test. 
One such case is the migration of the toroidal field at 
high latitudes. In dynamo wave (OL and IF) models 
at high latitudes all field components migrate pole- 
wards, following the Par ker-Yoshimura rule. In con- 
trast, the toroidal field in BL models typically shows 
an equatorward migration, advected by the return 
flow of meridional circulation near the bottom of the 
convective zone. Thus, if a clear signature of mi- 
grating high-latitude toroidal fields were found, this 
could solve the problem in either way, depending on 
the direction of migration. It is often claimed (e.g. 



Harvey 1994) that high latitude ephemeral active re- 



gions show an equatorward migration. But a closer 
examination of such claims shows that they are based 
only on the fact that high latitude ephemeral regions 
as a whole tend to lie in the backwards extension 
of the low-latitude butterfly "wings" (just as polar 
faculae do), and not on a detailed study of migra- 
tion patterns among high-latitude regions only. On 
the other hand, pallcbaut fc Makarov (1992) claim 
that at least 50 % of polar faculae (well known for 
their poleward drift) correspond to dipoles with a 
preferential east-west orientation, thus forming part 
of the toroidal field. It has even been claimed that 
the highest-latitude part of the sun spot butterfly di- 
agram also shows a poleward drift ( Becker 1959] ). A 



4.2. High-latitude field patterns 

In order to be able to decide between various dynamo 



clarification of this issue would clearly be important. 

A related problem is the origin of ephemeral ac- 
tive regions in general. If they are used as magnetic 
tracers we would obviously like to know what kind 
of field do they trace? 

4.3. Origin and profile of a 

It is sometimes claimed that the origin of the a- 
effect assumed is a basic difference between vari- 
ous dynamo models, so that excluding a certain a- 
mechanism amounts to excluding a dynamo variety. 
In reality, our knowledge of all possible types of a- 
effects is so scarce that wildly differing a-profiles can 
be derived for the same mechanism, depending on 
particular modelling assumptions. And vice versa: a 
given profile, used as input in a dynamo model, may 
be the consequence of various a-mechanisms. 

More in-depth study of the possible a-mechanisms 
that can reduce this uncertainty could nevertheless 
really be used to constrain dynamo models in the 
future. For now, I will argue that, independently of 
its mechanism, a should be expected to be concen- 
trated towards the bottom of the convective zone, or 
at least to pervade it uniformly (and not to be con- 
centrated to the surface, as assumed in BL models). 
Let us see this for the four a-mechanisms that have 
been proposed. 

- Cyclonic convection. In this classic variety 
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of the a-effect, proposed by Parker (|1955| ), a 
results from the effect of rotation on convec- 
tion, so it increases with the Coriolis number, 
peaking near the bottom of the solar convective 
zone. Its sign is positive in the unstable layer 
and negative in the overshoot layer below. 

Magnetostrophi c waves . Thi s mechanism 
was proposed by [Bchmitt ( 1987 ) specifically 
to suggest an alternative a-mechanism for the 
overshoot layer that will not be suppressed by 
the strong toroidal fields there. It consists 
in growing helical waves and its sign is neg- 
ative, as required by the dynamo wave models. 
More detailed calculations of its profile were 
performed by Ferriz Mas & Schiissler (1994), 



though the results may rely too heavily on 
the background stratification of the (incorrect) 
overshoot model used. 

Flux loop alpha. The magnetix flux loops 
emerging through the convective zone and cre- 
ating the active regions are essentially large- 
amplitude nonlinear versions of the above men- 
tioned unstable magnetostrophic waves, giving 
rise to a similar (but positive) a-effect. It 
is clear that after their decay, owing to their 
tilt, active regions contribute to the passive 
poloidal field. But is this contribution lim- 
ited to the surface, or do the "feet" of the 
loop similarly decay, down to the bottom of 
the convective zone? This problem relates to 
the next subsection, but one expects that the 
eroding action of external turbulence (Petrovay 
& Moreno-Insertis 1997) will lead to the de- 
cay of the whole loop, thereby extending the 
a-effect down to its footpoints. In fact, our 
turbulent erosion calculations show that 90% 
of the magnetic flux in an emerging flux loop 
is lost before it reaches the surface! This lost 
flux, already submitted to the action of Corio- 
lis force, should contribute to the poloidal field 
and to the a-effect far more than the actual 
active regions. And then we did not even men- 
tion the possibility of "failed active regions": 
flux loops that, their field strengths being too 



weak, never make it to the surface ( Petrovay fe 
Szakaly 1993 , Moreno-Insertis et al. 1995| ) 



All 

this strongly speaks for a flux loop a concen- 
trated to the bottom of the convective zone. 



Unstable Rossby waves. Recently, Dik 



pati & Gilman (2000) considered a shallow- 
water model of the tachocline, analysing 
the behaviour of small perturbations of its 
geostrophic equilibrium state. They find that 
such large-scale perturbations are unstable if 
the subadiabaticity is low enough and they are 



characterized by a correlation between the ver- 
tical components of velocity and vorticity, i.e. 
by a non-vanishing mean helicity. In other 
words, these perturbations essentially behave 
like global-scale unstable Rossby waves, and 
owing to their mean helicity they can be ex- 
pected to yield an a-effect. 

4.4. The ultimate fate of emerged flux 

This problem is interesting in its own right, as well as 
because of its importance for the a-effect problem. 
Once the active region has decayed, what happens 
to the "trunks" of the magnetic trees? Will they 
just stay there, "bleeding"? Will they reconnect and 
be drawn back below the convective zone? Or will 
they also decay? As indicated above, the turbulent 



erosion models ( Petrovay fc Moreno- Inserts 1997 ) 
seem to support the latter possibility. At any rate, 
thorough observational studies of the active region 
decay (cf. the review by van Driel-Gesztelyi 1999| ) 
may shed more light on this problem. 



5. CONCLUSION: ...AND THE 
SIMULATIONS? 



As we have seen, there is still no really convincing 
answer to the question "What makes the Sun tick?" . 
Several conflicting approaches exist, none of which 
can fully explain the observed features of the solar 
cycle and all of which rely on some more or less arbi- 
trary assumptions. In contrast to the optimistic out- 
look of about three decades ago, nothing seems safe 
now, not even the dynamo wave origin of the sunspot 
butterfly diagram. And the prospects do not seem 
to promise a spectacular change in this situation in 
the near future. 



A number of previous reviews (e.g. Weiss 1994) 
expressed the view that the ultimate solution to the 
solar dynamo problem should be expected from nu- 
merical hydrodynamical simulations. So the reader 
may ask: why was the simulational approach hardly 
even mentioned in this review? Why is it that only 
mean field models were treated? The answer does 
not lie in the (undeniable) bias of the author but in 
the fact that in the past decade there simply have not 
been any numerical simulations constructed with the 
ambition of producing a global model for the solar 
cycle. As we indicated above, such attempts were 
made back in the 1980's; but, after their failure in 
reproducing the butterfly diagram and the internal 
rotation profile became clear, this line of research 
was completely abandoned. Those very few hydro- 
dynamical simulations in the 1990's that had some 
direct bearing on the solar dynamo problem aimed at 
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the modelling of smaller volumes inside the connec- 
tive zone with the main purpose of studying the fine 



str ucture and tr a nsport of the magnetic field ( Nord 



lun i et al. 1992 , Tobias et al 
lun|d| E000|) . 



1998, Dorch & Nord- 



The belief in the "miraculous healing power" of 
simulations may indeed have been too zealous be- 
fore. After all, even in an ideal future where an in- 
finite computing potential would make it feasible to 
make a 3D direct numerical simulation of the whole 
Sun for 4.6 • 10 9 years, the result would have no more 
direct benefit for us (apart from demonstrating that 
the laws of classical magnetohydrodynamics are in- 
deed sufficient to describe solar phenomena) than of- 
fering a chance to determine any physical quantity 
at any internal point of our star with arbitrary pre- 
cision. This feat, though far beyond the possibilities 
of contemporary observational solar physics, is still 
basically experimental work that, though important, 
in no way can substitute "real" theory. 

And yet, the other extreme, represented by the 
present complete abandonment of the simulational 
approach, is just as lamentable as its opposite. This 
is especially so as the great advance in the geody- 
namo problem I referred to in the Introduction can 
mainly be attributed to the success of numerical sim- 
ulations. The following pseudo-historical quote by 
Borges fits nicely the present situation in this re- 
spect. 



Jorge Luis Borges: 

On Rigour in Science 

...In that Empire, the Art of Cartography achieved 
such Perfection that the map of a single Province oc- 
cupied a whole City, and the map of the Empire, a 
whole Province. With time, these Gigantic Maps 
proved unsatisfactory and the Colleges of Cartog- 
raphers set up a Map of the Empire that had the 
size of the Empire and coincided with it point by 
point. Less addicted to the Study of Cartography, 
the Following Generations considered this extensive 
Map useless and, not without Disrespect, they aban- 
doned it to the Mercy of the Sun and of the Seasons. 
In the Western deserts broken Ruins of the Map still 
persist, inhabited by Animals and Beggars; in all the 
Country no other Relicts of the Geographic Disci- 
plines remain. 

Suarez Miranda: Viajes de varones prudentes, 
libro cuarto, cap.XLV, Lerida, 1658 
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